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We study the quantum phase transitions in the two-dimensional spin-orbit models in terms of 
fidelity susceptibility and reduced fidelity susceptibility. An order-to-order phase transition is iden- 
tified by fidelity susceptibility in the two-dimensional Heisenberg XXZ model with Dzyaloshinsky- 
Moriya interaction on a square lattice. The finite size scaling of fidelity susceptibility shows a 
power-law divergence at criticality, which indicates the quantum phase transition is of second order. 
Two distinct types of quantum phase transitions are witnessed by fidelity susceptibility in Kitaev- 
Heisenberg model on a hexagonal lattice. We exploit the symmetry of two-dimensional quantum 
compass model, and obtain a simple analytic expression of reduced fidelity susceptibility. Com- 
pared with the derivative of ground-state energy, the fidelity susceptibility is a bit more sensitive 
to phase transition. The violation of power-law behavior for the scaling of reduced fidelity suscep- 
tibility at criticality suggests that the quantum phase transition belongs to a first-order transition. 
We conclude that fidelity susceptibility and reduced fidelity susceptibility show great advantage to 
characterize diverse quantum phase transitions in spin-orbit models. 

PACS numbers: 03.67.-a,64.70.Tg,75.25.Dk,05.70. Jk 



I. INTRODUCTION 

Spin-spin interactions have been intensively studied in 
quantum magnets and Mott insulators in the last decades 
[T]. The effect of the orbital degree of freedom has re- 
ceived much attention since the discovery of a variety of 
novel physical phenomena and a diversity of new phases 
in transition metal oxides (TMOs) [2l[3j. In particular, 
under an octahedral environment, the d orbital degen- 
eracy of transition metal ions is partially lifted, and the 
remaining orbital degrees of freedom can be generally 
described by localized S = 1/2 pseudospins. A typical 
effect induced by such a symmetry degradation of or- 
bital degeneracy is the presence of bond-selective pseu- 
dospin interaction. It is because the spatial orientations 
of the orbits lead to anisotropic overlaps between neigh- 
boring ions. Consequently, the interactions among differ- 
ent bonds are intrinsically frustrated. To understand the 
orbital degree of freedom, the so-called quantum compass 
model (QCM) has been given rise to intensive research 
j4HT5]. In QCM, the pseudospin operators are coupled in 
such a way as to mimic the built-in competition between 
the orbital orders in different directions [H |3] . The frus- 
tration leads to macroscopic degeneracy in the classical 
ground state [4, and even highly degenerate quantum 
ground state [131 E] • Interestingly, the two-dimensional 
(2D) QCM has become a prototype to generate topolog- 
ically protected qubits [TS] |T6] . 

Besides, considering the spin-orbit coupling in TMOs 
is inevitable and intriguing. For example, octahedra tilt 
may give rise to effective Dzyaloshinsky-Moriya interac- 
tion (DMI) [17]. Comparing those that only possess spin- 
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spin coupling, spin-orbit models appear to be more intri- 
cate. The interplays between spin and orbital degree of 
freedom host a variety of different phases. For instance, 
the orbital exchange in a honeycomb lattice induces or- 
bital ordering [18l[T9] and topological order [20 . 

In this paper, we concentrate on portraying the quan- 
tum phase transitions (QPTs) in 2D spin-orbit interac- 
tion Hamiltonians. As we know, a QPT identifies any 
point of nonanalyticity in the ground-state (GS) energy 
of an infinite lattice system [21 . Conventionally, local 
order parameters are needed to detect the nonanalytic- 
ity in the GS properties as the system varies across the 
quantum critical point (QCP). However, the knowledge 
of the local order parameter is not easy to retrieve from 
a general many-body system, especially for QPTs be- 
yond the framework of the Landau- Ginzburg symmetry 
breaking paradigm [20l [22] . Recently, quantum fidelity, 
also referred to as the GS fidelity, sparked great interest 
among the community to use it as a probe for the QCP 
[23l [24] , The fidelity defines the overlap between two 
neighboring ground states of a quantum Hamiltonian in 
the parameter space, i.e., 

F(A,(5A) = K^o(A)|^o(A + a))|, (1) 

where |^o(A)) is the GS wave function of a many-body 
Hamiltonian H{\)=Hq + Ai^/, A is the external driv- 
ing parameter, and 5\ is a tiny variation of the external 
parameter. Though borrowed from the quantum infor- 
mation theory, fidelity has been proved to be a useful 
and powerful tool to detect and characterize QPTs in 
condensed matter physics [25 . In order to remove the 
artificial variation of external parameters, one of the au- 
thors and collaborators in Ref. [26 introduced the con- 
cept of fidelity susceptibility (FS), which is the leading- 
order term of fidelity with respect to the external driving 
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parameter, 
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FS elucidates the rate of change of fidehty under an in- 
finitesimal variation of the driving parameter. There ex- 
ists an intrinsic relation between the FS and the deriva- 
tives of GS energy as following: 
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Here the eigenstates |^n(A)) satisfy i^(A)|^n(A)) = 
En{X)\'^nM)' The apparent similarity between Eq. (|3| 
and Eq. Q arouses similar critical behavior around the 
critical point. 

In addition, the reduced fidelity and its susceptibil- 
ity were also suggested in the studies of QPTs [27H3Q] . 
The reduced fidelity concerns the similarity of a local re- 
gion of the system with respect to the driving parameter. 
Despite the locality, the reduce fidelity fidelity (RFS) en- 
codes the fingerprint of QPTs of the whole system, and is 
even more sensitive than its global counterparts [311 [32] . 

For a second-order QPT, around the critical point, the 
correlation length ^ diverges as (A — Ac)~^, while the gap 
in the excitation spectrum vanishes as (A — Ac)^^, where 
Ac is the critical point in the thermodynamic limit, u 
and z are the critical exponents. At the critical point, 
i.e., ^ = oo, the only length scale is the system size L. 
We can account for the divergence at the QPT by formu- 
lating a finite-size scaling (ESS) theory. Universal infor- 
mation could be decoded from the scaling behavior of FS 
[33fl35j . FS increases as the system size grows, and the 
summation in Eq.([3| contributes to an extensive scaling 
of Xf in the off-critical region. Therefore, FS per site 
Xf/N appears to be a well-defined value, where N = 
is the number of sites and d stands for system dimen- 
sionality. Instead, FS exhibits stronger dependence on 
system size across the critical point than in non-critical 
region, showing that a singularity emerges in the sum- 
mation of Eq.([3|. This implies an abrupt change in the 
ground state of the system at the QCP in the thermo- 
dynamic limit. Following standard arguments in scaling 
analysis [36], one obtains that FS per site scales as ^37i - 
[39] 



Xf/N - L^/^ 
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Due to the arbitrariness of relevance of the driving Hamil- 
tonian Hj under the renormalization group transforma- 
tion, Xf/N could be (i) superextensive livd <2 [37], (ii) 
extensive \i vd = 2 [JO] [41] , (iii) subextensive if ud > 2 
[42]. From Eq. ^ we can deduce that the ground state 
of the system should be gapless if xf is superextensive, 
but not vise versa. For finite size system, the position 



of a divergence peak defines a pseudocritical point X^ as 
the precursor of a QPT, and it approaches the critical 
point Ac as L ^ oo. For sufficiently relevant perturba- 
tions on large size system, i.e., ud < 2, the leading term 
in expansion of pseudocritical point obeys such scaling 
behavior as [43l 
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Thus, the behavior of xf on finite systems in the vicinity 
of a second-order QCP can be estimated as 
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where / is an unknown regular scaling function, and C 
is a constant independent of A and L. 

However, many important QPTs fall in the cate- 
gory of first-order QPTs, in which the first derivative 
of GS energy exhibits discontinuity at critical point. 
Different from second-order QPT, there is no charac- 
teristic correlation length ^ in first-order QPT, and 
in general the FSS will violate the scaling relations 
of second-order QPTs [47]. For a typical first-order 
QPT, two competing ground states |^<(A)) and |^>(A)) 
are degenerate at critical point Ac in the thermody- 
namic limit, and they become energetically favorable 
on one side of Ac such as |^<(A)) = \^^(X < Ac)) 
and |^>(A)) = |^o(A > Ac)) respectively. The level 
crossings in the thermodynamic limit usually turn into 
avoided level crossings for finite system, and degener- 
acy at critical point in general is lifted, opening an en- 
ergy gap A^. In the low-energy subspaces spanned by 
|^<(Ac)) and |^>(Ac)), the diagonal matrix elements 
coincide (^<(Ac)|i/|^<(Ac)) = (^>(Ac)|i/|^>(Ac)), but 
the off-diagonal matrix elements (^<(Ac)|i^|^>(Ac)) = 
((^>(Ac)|i^|^<(Ac)))* induce an avoided level crossing 
with 



Ag = 2|(vI/<(Ae)|if|M/>(Ae))|. 



(8) 



For a Hamiltonian with local interactions (e.g., nearest 
neighbors). 
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|^>(Ac)) = T^nes ^ril^n) , \^m) (|^n)) form a complete 
set of basis vectors in Hilbert space 5, and Cm (cn) are the 
corresponding coefficients. The dimension of S is of 
{V is the degree of freedom of Hamiltonian constituent, 
e.g., P = 2 for spin- 1/2 Hamiltonian), and Cm (cn) are of 
order The finite Hamming distance between \(pm) 

and \(pn) gives limited number (roughly speaking, 0{N)) 
of nonzero H^n. Hence, the off-diagonal matrix elements 
scale exponentially (^<(Ac)|i^|^>(Ac)) - 0{NV-^), 
and consequently - 0{NV-^) [48] [49]. If any of the 
ground states are degenerate protected by symmetries, 
H should be written in the subspace of all the degener- 
ate states, and the dimension of the matrix will be larger 
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than 2. Considering Eq. ([3|, the exponentiahy closing 
gap gives a dominant contribution, manifesting that the 
FS xf should carry the signature of an exponential di- 
vergence at critical point. Moreover, from another point 
of view, due to their macroscopic distinguishability, the 
overlap between states |(^<(A)|^>(A))| should be expo- 
nentially small with system size N at criticality [48 . In 
general, FS per site scales exponentially, i.e.. 



Xf/N - i/(iV)e-^^, 



(10) 



where /i is a size independent constant, and g{N) is a 
polynomial function of which is a correction to the 
exponential term. 

In this work, we use FS and RFS to incarnate the crit- 
ical phenomena in 2D spin-orbit models. We show that 
the FS and RFS manifest themselves as extreme points 
at the QCPs. In other words, FS and RFS are quite sen- 
sitive detectors of the QPTs in spin-orbit models. This 
article is organized as follows: In Sec. [llj a general Hamil- 
tonian is given. We investigate QPTs in the spin-orbit 
model on a square lattice in Sec. |III[ The phase dia- 



gram and the FSS are studied through FS. Next, in Sec. 
JVj we study QPTs in Kitaev-Heisenberg model on a 2D 
honeycomb lattice by FS and the second derivative of the 
the GS energy. Both approaches identify three distinct 
phases, and FS behaves more sensitively. After that, we 
take advantage of the RFS to locate the QCP of 2D com- 
pass model in Sec. [V| and show that the exponentially 
divergent peaks of FS and RFS imply that the first order 
QPT. A brief summary is presented in Sec. VI 



II. GENERAL HAMILTONIAN 

Motivated by a broad interest in Mott insulators, we 
construct the general structure of model Hamiltonian on 
a given 7 bond connecting two nearest-neighbor (NN) 
sites i and j. 



^ij - ^7^i 
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The first term in Hamiltonian (11) describes bond- 
selective interaction, where the orbital exchange con- 
stant J J = At^/U is derived from multi-orbital Hub- 
bard Hamiltonian consisting of the local on-site repul- 
sion U and the hopping integral t^. The second term 
corresponds to the isotropic Heisenberg coupling, and 
Dzyaloshinsky-Moriya anisotropy D in the third term 
comes from the lattice distortions. The Hamiltonian ( 11 ) 



extrapolates from the Heisenberg model to QCM depend- 
ing on the bonds geometry. In this paper, we devote our 
study to 2D lattices, either square or hexagonal lattice. 
Several real systems motivate the investigation of these 
kinds of lattices, such as SrTiOs gO], Sr2lr04 [5TH53]. 
The flexibility of parameters induces a rich variety of 
phases and various fascinating physical phenomena. 



III. FIDELITY SUSCEPTIBILITY IN THE 
SPIN-ORBIT MODEL 

In strong spin-orbit coupling materials, an XXZ model 
with DMI was suggested. 
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where A is the anisotropic parameter and D is the 
strength of DMI. A general boundary condition can be 
written as S^j^-^ = pS"^ ^ where = S? ±iSy , p = cor- 
responds to open boundary condition (OBC) and p = 1 
to the periodic boundary condition (PBC). The Hamil- 
tonian (12) was proposed to describe the layered com- 
pound Sr2lr04 [17]. The spin canting is induced by 
lattice distortion of corner-shared IrOe octahedra. The 
DMI was introduced originally to explain the presence of 
weak ferromagnetism in antiferromagnetic (AFM) mate- 
rials ^541 [55], such as a-Fe2 03, MnCOs and C0CO3, since 
such antisymmetric interaction could produce small spin 
cantings. Recently, the influence of DMI has become very 
important in elucidating many interesting properties of 
different systems, e.g., ferroelectric polarization in multi- 
ferroic materials [56: , exchange bias effects in perovskites 
[57] . asymmetric spin- wave dispersion in double layer Fe 
[58] [59] , and noncollinear magnetism in FePt alloy films 
[60] , For simplicity, we assume the DMI fluctuation is 
acting on the x-y plane, and the vector D is imposed 
along the z direction, i.e., D = Dz. 

The Hamiltonian (12) reduces to the anisotropic 
Heisenberg XXZ model when the rotations of IrO^ oc- 
tahedra are absent, i.e., D=0. When A ^ 1, one- 
dimensional (ID) spin- 1/2 XXZ chain has long range or- 
der and gapped domain-wall excitations. On the other 
hand, in the XX limit, i.e., A = 0, the system is equiva- 
lent to a chain of non-interacting fermionic model, which 
becomes gapless in the thermodynamic limit. The QPT 
takes place at the isotropic point Ac = 1, which is a 
Berezinskii-Kosterlitz-Thouless (BKT) transition point. 
BKT phase transition belongs to an infinite-order phase 
transition, and the ground-state energy and all of its 
derivatives with respect to A are continuous at the crit- 
ical point. However, the FS succeeds in detecting the 
nonanalyticities of the ground state across BKT transi- 
tion [61]. For instance, the BKT transition from gapless 
Tomonaga-Luttinger liquid to gapped Ising phase in ID 
XXZ model is detected by the divergence of the FS us- 
ing the density-matrix-renormalization-group (DMRG) 
technique [62 . In addition, the BKT transition from 
spin fluid to dimerized phase in the J1-J2 model [63 and 
the superfiuid-insulator transition in the Bose-Hubbard 
model at integer filling [35 are also able to be signaled 
by FS. 

Different from the ID case, 2D XXZ model exhibits 
a second-order QPT at the isotropic point Ac=l in the 
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thermodynamic limit, where the first excited energy lev- 
els cross [64 . For A 1, the Ising term in the Hamil- 
tonian dominates and the ground state is an AFM phase 
along the z direction. For A <C 1, the first two terms in 
the Hamiltonian dominate and the ground state is also 
an AFM phase, but in the x-y plane. It is well known 
that long-range orders are present in both phases [65\ 
FS can serve as a sensitive detector of the critical point 
in the 2D XXZ model on a relatively small square lattice 
[34 . 

Actually, the DMI does not change the universality 
class of XXZ model. The DMI can be eliminated from 



the Hamiltonian (12) by a spin axes rotation [66r;68 . and 
after rotation an unitarily equivalent form is given by 



Hdm(A,I)) 



cos ( 
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A — Acos(/) and Hxxz(A) = 
-AS^S? with boundary condition 

Note that the mapping becomes 
exactly equivalent only in the thermodynamic limit and 
for open boundary condition. In the thermodynamic 
limit, the boundary condition does not affect the critical 
behavior and consequently the 1-Lbm{^^ D) will have the 
same critical properties as the Hxxz(A)- Hence, QCP 
at Ac=l in Hxxz becomes a critical line Ac = Vl + 1)^ 
in Hdm. 

FS and the second derivative of GS energy of N = 20 
square lattice (see Fig. [T]) as a function of anisotropy A 
for D = and D = 2 are demonstrated in Fig. |2j As 
is sketched, both approaches display broad peaks at an 
identical A. However, the peak of FS is slightly narrower 
than that of second derivative of GS energy. The loca- 
tions of peaks are specified as pseudocritical points. The 
heights of local maxima decrease as D increases. 

Next, to check whether these extreme points can be re- 
garded as QCPs, we should perform FSS analysis, as is il- 
lustrated in Fig. |3] In Figs, [sj^a) andjsjb), we plot FS per 
site as a function of A for different lattice size A/'=10, 16, 
18, 20 and 26. As is clearly shown that xf becomes more 
pronounced for increasing N. Xf is extensive in the off- 
critical region, while superextensive at the pseudocritical 
point. From the scaling relation Eq. ([5|, a linear depen- 
dence for the maximum fidelity susceptibilities Xf^^ on 
L^/^ is expected with effective length L = \/N . This is 
confirmed by the results shown in Figs, [sjc) andjsjd), in 
which we plot the logarithm of maximum fidelity suscep- 
tibilities Xf pseudocritical points verse a function of 
hiN. For large Inxp^^ scales linearly with In A/", and 
the slope is interpreted as the inverse of critical exponent 
of the correlation length l/u from Eq. ([5|. By apply- 
ing linear regression to the raw data obtained from XXZ 
model on various square lattices, we derive 1/v ^ 2.68(6) 
for = 0, while l/v ^ 2.62(8) for D = 2. The values 
got from the measurements oi D = {) and D = 2 are con- 
sistent with each other up to two digits. The validation 
of Eq. ^ implies that the QPT in 2D XXZ with DMI 




FIG. 1: (Color online) Two-dimensional square structures for 
different system sizes N=10, 16, 18,20,26. 



manifests itself as a clear sign of a second-order QPT. 

From Figs. |3|a) and [sj^b), the positions of cusps 
seemingly converge toward the critical points. How- 
ever, these system sizes seem insufficiently to estimate 
the critical exponents from Eqs. (|6| and ([7|), and the 
corrections to scaling relations are not negligible. In 
this case, the shift of the location of the pseudocriti- 
cal point from real critical point should be replaced with 
\Al-Ac\ ^ ciL"^/^ + C2L"^/^ + . . ., where ci and C2 are 
the coefficients [43 . As is shown in insets of Figs, [sj^c) 
andjsjd), we notice that the data points corresponding 
to the small system sizes clearly deviate from the lin- 
ear fit obtained for the points for the two largest N. As 
A/" ^ oo, the pseudocritical points AJ^ approach the crit- 
ical points Ac- The linear fits from results of N = 20 
and = 26 (see Fig. [T]) yield the estimates the Ac are 
1.02 and 2.47, respectively. To get more precise critical 
exponents, more elaborate schemes are needed, such as 
large-scale quantum Monte Carlo simulations [371 HI] • 

To proceed, we calculate the FS on N=20 square lat- 
tice for different sets of parameters A and D. The FS as 
a function of I) for A = 1.5 and A = 2.0 are depicted in 



FIG. 2: (Color online) The second derivative of ground-state 
energy (solid line) and fidelity susceptibility (dash line) of 
N = 20 square lattice as a function of A for (a) D = and 
(b) D = 2. 



Fig. [4) and the FS as a function of A for D=0, 1,2, and 
3 are plotted in Fig. [5] A boost of A and D suppress the 
FS. By retrieving the critical points from the locations of 
FS peaks, we are able to draw the corresponding phase 
diagrams. The pseudocritical lines in the {D, A) plane 
and the (A, D) plane are shown in the insets of Fig. |4] 
and Fig. |5j We can notice that pseudocritical line for 
= 20 (solid dotted line) are qualitatively similar to 
the critical line in the thermodynamic limit (dash line). 



IV. FIDELITY SUSCEPTIBILITY IN THE 2D 
KITAEV-HEISENBERG MODEL 

Since Kitaev introduced a spin 1/2 quantum lattice 
model with Abelian and non-Abelian topological phases, 
finding a physical realization of Kitaev model has trig- 
gered a tremendous amount of interest [69]. It is pro- 
posed that in iridium oxides A2lr03, the strong spin-orbit 
coupling may lead to the desired anisotropy of the spin in- 
teraction. The Ir^+ ions in iridium oxides A2lr03 can be 
effectively illustrated as spin half on a honeycomb lattice, 
where three distinct types of NN bonds are referred to as 



FIG. 3: (Color online) Top: the fidelity susceptibility per site 
Xf/A as a function of anisotropy parameter A for various size 
N= 10, 16, 18, 20, 26 with (a) D = and (b) D = 2. Bottom: 
finite size scaling analysis for the height and location of the 
peaks. The logarithm of maximum fidelity susceptibilities 
Xf^^ for (c) Z)=0 and (d) D—2^ respectively, are plotted as a 
function of In A. The dash lines are least square straight line 
fits with l/v ^ 2.68(6) for (c) D=0 and l/u ^ 2.62(8) for (d) 
D—2. The insets in panels (c) and (d) show pseudocritical 
points A|^ for various sizes A correspondingly. 



7 = (0°, 120°, 240°) bonds. To describe the competition 
between direct exchange and superexchange mechanisms, 
a so-called Kitaev-Heisenberg model was dedicated, i.e.. 



-2aS7S] 



(l-a)Si-Sj. 



(14) 



<ij>ll7 



The Hamiltonian ( 14 ) has been parameterized, which en- 



compasses a few well-known models. For a = 0, it re- 
duces to the Heisenberg model on a hexagonal lattice, 
while it becomes exactly solvable ferromagnetic Kitaev 
model at a = 1. Another solvable point corresponds to 
a = 1/2, where it is unitarily equivalent to ferromagnetic 
Heisenberg model. The ground state of Hamiltonian ( 14) 
evolves from 2D Neel AFM state {a = 0) to stripy an- 
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FIG. 4: (Color online) The fidelity susceptibility xf as a 
function of D for A = 1.5 and A = 2.0. Inset shows the 
pseudocritical line retrieved from the peaks of fidelity suscep- 
tibility for N — 20 (solid dotted line) and the critical line 
Dc = VA^ — 1 in thermodynamic limit (dash line). 

tiferromagnetism (a = 1/2), and to spin liquid (a = 1) 
as a increases. We perform exact diagonalization (ED) 
on two geometries with A^=16 and A/'=24 (see Fig. [6| to 
calculate the GS energy Eq and fidelity susceptibility xf- 

The second derivative of energy density eo = —Eq/N 
and FS per site xf/^ as a function of coupling strength a 
are obtained by Lanczos calculation. As is illustrated in 
Fig. [7| two anomalies emerge when a increases from to 
1, indicating the system undergoes two phase transitions. 
The characterization of QPTs by FS is compatible with 
the second derivative of energy density. We find that they 
play similar roles in identifying the QPTs |63]. However, 
the former demonstrates more pronounced peaks than 
the latter (notice the y-axis is of a logarithmic scale). 
QPT from Neel-ordered to stripe-ordered phase takes 
place at a ^ 0.4, and this order-to-order phase transition 
seems insensitive to system size N for these two configu- 
rations. On the other hand, the second order-to-disorder 
phase transition at a ^ 0.8 appears to be sensitive to 
N. The numerical calculations point out that the QPT 
around a = 0.4 is first-order, while there is no consensus 
on the character of the QPT around a = 0.8 £0l 171] . 




FIG. 5: (Color online) The fidelity susceptibility xf in the 
ground state of the 2D 20-site XXZ model in terms of 
anisotropy for different values of DMI. Inset shows the pseudo 
critical line by the extraction from the maximum points of the 
fidelity susceptibility (solid dotted line), and the critical line 
Ac = a/1 + in thermodynamic limit (dash line). 

However, it is difficult to carry on ED on larger cluster 
for the hexagonal geometry and substantiate the scaling 
behavior of 2D Kitaev-Heisenberg model unless sophisti- 
cated techniques are used ^37l [44] . 



V. REDUCED FIDELITY SUSCEPTIBILITY IN 
THE 2D COMPASS MODEL 

In recent years, the 2D AFM QCM has attracted con- 
siderable attention due to its interdisciplinary character. 
On one hand, it plays an important role in describing 
orbital interactions in TMOs. On the other hand, 2D 
QCM is equivalent to Xu-Moore model [72j and toric 
code model in a transverse field ^73^, which could pos- 
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16 sites 



24 sites 




FIG. 6: (Color online) Two-dimensional structures for system 
sizes N—16 and 24, which can be placed on a hexagonal lattice 
with periodic boundary conditions. 



sibly be used to generate protected qubits realized by 
Josephson-coupled p±ip superconducting arrays [15j . 

QCM is defined on a = L x L square lattice with 
PBC by the Hamiltonian 



(15) 



(iJ>ll^ 



(iJ>l|ez 



where Gx, are unit vectors along the x and z directions, 
and Jx (Jz) is the coupling in the x (z) direction. 
(a=x, z) are the pseudospin operators of lattice site i 
obeying [S^^ S?] = icafB-^Sy Si^y There exists a first-order 
phase transition at the self-dual point, i.e., Jx = Jz^ and 
it is widely believed that the phase transition belongs to 
first-order [Ml [73]. 

In the case of L being even, this model is equivalent 
to the ferromagnetic QCM by rotating the pseudospin 
operators at one sublattice by an angle tt about iS^-axis. 
The Hamiltonian ( 15 ) enjoys such commutation relation 
with the parity operators Pj and Q^, which are defined 
as 
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n (2^m) = n ^iv Q^=i[ i^sf,) = n 



(16) 

where indices i and j are the x- and z-component of lat- 
tice site i. However, column parity operator Pj does 
not commute but anticommutes with row parity oper- 
ator Qi. In such circumstance, the Hilbert space can 
then be decomposed into subspace V{{pj})^ where pj is 
the eigenvalue of Pj and is specified in each subspace. 
Since Pj = 1, pj can be either 1 or -1. The lowest energy 
of the Hamiltonian in the subspace V{{pj}) is nonde- 
generate. In Ref. [74], the authors unraveled the hidden 
dimer order therein. The symmetries in these systems in- 
duce large degeneracies in their energy spectra, and make 
their numerical simulation tricky. Recently, the ground 
state of Hqcm is proven to reside in the most homoge- 
neous subspaces V{{pj = 1}) and V{{pj = —1}), in other 
words, they are two- fold degenerate [75 . 

At this stage, a single-site pseudospin- flipping interac- 
tion will change the parity of the chain, contrary to the 
flipping terms in H. For example, or acting on 
one chain along x direction will change the parity of the 
chain. For two-point correlation functions between sites 
i and j separated by n sites, it is not difficult to yield 



(-1)" vl/o(A) S?S?+„, *o(A)) = 0. (17) 



FIG. 7: (Color online) (a) The second derivative of energy 
density eo verse a for N = 16 and N = 24. (b) Fidelity 
susceptibility per site Xf/N as a function of a. 



Similarly, e^^=0, e^^ = 0, e^^ = 0. The only correlation 
functions surviving are Cnx = 4e^^ and Cnz = ^^nz^ 
are shown in Fig. [Sjand inset of Fig. 11 As Jx increases, 
Cz decreases while Cx increases accordingly. They become 
equal at Jx= Jz- 

With these properties of the ground state |^o) of the 
2D QCM, a simple form for the two-site density matrix 
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is obtained, 



p{i,j)=Tr'{\^o){^o\) = \ ^ (<tX')<t,X', (18) 

a,a'=0 

in which the prime means tracing over ah the other pseu- 
dospin degrees of freedom except the two sites i and j. 
cr" are PauU matrices a^, and for a = 1 to 3, and 
2 by 2 unit matrix for a=0. If the two pseudospins are 
Unked by an x-type bond, we have 

p,(i,j) = Tr'(|*o>{*o|) = HS?S?)S?S? + i/i/j, (19) 

in which Ij and /j are the 2 by 2 unit matrices. If trans- 
lational invariance is preserved in the ground state l^o), 
the reduced density matrix can be simphfied as 



(20) 



Similarly, if the two sites i and j are linked by z-type 
bond, we have 



p,(i,j) = -C,S?S? + -IJy 



(21) 



In addition, the NN two-point correlation functions Cq 
can be calculated using the Feymann-Hellmann Theorem 



4 dEo 



(22) 



where a = x^z means x-type and z-type bond respec- 
tively. 

The reduced fidelity Fr is defined as the overlap be- 
tween paihi) and Paihi), i-e. 



1/2 , 1/2 



(23) 



The prime means a different reduced density matrix in- 
duced by the tiny changes of driving parameters in the 
Hamiltonian. By Eqs. (20) and (21), it is obvious that 
pa commutes with p'^, and we can easily evaluate the 
reduced fidelity 

FMhj),p'a(hj)) = ^(V(l + Ca)(l+C^,) 

+ v/(l-C)(l-C^)). (24) 

The two-site RFS Xr could also be obtained straightfor- 
wardly. 



-2\nFr _ {dj^Ca? 



(25) 



in which Jjs [/3=x,z) is the driving parameter for the 
QPT. According to Eqs. (p2|) and ([25|, we observe that 



the numerator of Xr"^ is proportional to the square of the 
second derivative of GS energy and the denominator is 




FIG. 8: (Color online) The nearest-neighbor correlation func- 
tions Cx and Cz with respect to Jx/Jz for = 16 square 
lattice. Two curves cross at Jx / J z — 1 • 



finite in unpolarized state. The second power in the nu- 
merator indicates that the two-site RFS is more effective 
than the second derivative of the GS energy in measuring 
QPTs [76 . We compare the RFS of two NN sites with 
the second derivative of GS energy of = 18 square 
lattice (see Fig. [T]) displayed in Fig. [9j and find that 
they indeed present cusp-shaped peaks, but the former 
is more sensitive to phase transition than the latter. The 
pronounced maxima of peaks arise at Jx/Jz = 1- 
Fig. 
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unveils the results from ED on various 2D 
square lattices with PBC. As the coupling J^ changes, 
both FS and RFS of two NN sites exhibit peaks around 
Jx = Jz- With increasing system size, the peaks of xf (or 
Xr) become more pronounced and pseudocritical points 
seemingly converge toward the real critical point Jx = Jz 
quickly. We plot the maximum fidelity Xf maxi- 
mum reduced fidelity susceptibility x^^^ against the sys- 
tem size N in Fig. [lOjd). The scaling of the peaks re- 
veals approximately exponential divergences at critical- 
ity instead of power-law divergences, which indicates the 
occurrence of non-second-order phase transition. 

As is argued in Eqs. ([8| - (10), the approximately 
exponential divergence of FS and RFS hints an exponen- 
tially small gap at the critical point and thus a first-order 
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FIG. 9: (Color online) The second derivative of ground state 
energy (solid line) and reduced fidelity susceptibility (dash 
line) of A/" = 18 square lattice as a function of Jx/ Jz- 



phase transition. For two next-nearest-neighbor (NNN) 
sites, we calculate the RFS according to Eq. (25). We 



observe a similar behavior of the RFS with respect to 
Jx/ Jzi and find that the RFS of NNN pair is a bit larger 
than the RFS of NN bond. This is illustrated in Fig. |TT 
In a word, the two-site RFS can serve as a signature for 
the QPTs in the 2D compass model. 



VI. SUMMARY AND DISCUSSION 

Motivated by recent theoretical and experimental work 
on orbital degree of freedom in Mott insulators, we have 
studied the QPTs in various 2D spin-orbit models us- 
ing FS and RFS. The numerical analysis is performed on 
2D clusters with the Lanczos algorithm. The spin-orbit 
model hosts plentiful phases, including phase transitions 
within and beyond the framework of Landau-Ginzburg 
paradigm. For the 2D XXZ model with DMI, with 
increasing driving parameters, the system undergoes a 
second-order phase transition from AFM state along x-y 
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FIG. 10: (Color online) (a) The fidelity susceptibility and (b) 
the reduced fidelity susceptibility of x-bond and (c) 2;-bond in 
the ground state of the 2D compass model on a square lattice 
as a function of Jx/Jz- The curves correspond to different 
lattice sizes A/'=8, 16 and 18, respectively, (d) The log-linear 
plot of maximum fidelity susceptibility and maximum reduced 
fidelity susceptibility against the system size N . 



plane to AFM state along z direction. We compare FS 
with second derivative of GS energy, and find both of 
them exhibit similar peaks. The FSS demonstrates that 
FS per site should diverge in the thermodynamic limit at 
pseudocritical point, and the locations of extreme points 
approach QCPs accordingly. The power-law divergence 
of FS at criticality indicates the quantum phase transi- 
tion is of second order and the critical exponent v is ob- 
tained. Analogously, as the exchange coupling changes, 
the ground state of 2D Kitaev-Heisenberg model evolves 
from Neel AFM state to stripy AFM state, and to Kitaev 
spin liquid. The QCPs could be signaled by the peaks of 
both FS and second derivative of GS energy. The non- 
local symmetries in 2D AFM QCM guarantee that the 
RFS can be written in an analytical form as a function 
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of the correlation functions. Peaks of FS and two-site 
RFS take place around Jx=Jz' The quasicritical point 
develops into critical point quickly with increasing the 
system size. Scaling of the peaks reveals an exponential 
divergence at criticality, which suggests that a first-order 
phase transition happens. In conclusion, the FS and RFS 
are effective tools in detecting diverse QPTs in 2D spin- 
orbit models, and their scaling behaviors may hint the 
orders of phase transitions. 
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FIG. 11: (Color online) The reduced fidelity susceptibility 
of two next-nearest-neighbor sites for = 18 square lattice 
with respect to Jx/ Jz- Inset shows the next-nearest-neighbor 
correlation functions verse Jx j J z- 
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